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^ ABSTRACT 

O [ Context. The first couple of stellar generations may have been massive, of order 100 M Q , and to have played a dominant role in 

D 1 galaxy formation and the chemical enrichment of the early Universe. Some fraction of these objects may have died as pair-instability 

, supernovae or gamma-ray bursts. The winds of these stars may have played an important role in determining these outcomes. As 

the winds are driven by radiation pressure on spectral lines, their strengths are expected to vary with metallicity. Until now, most 
mass-loss predictions for metal-poor O-type stars have assumed a scaled-down solar-abundance pattern. However, Population III 
evolutionary tracks show significant surface enrichment through rotational mixing of CNO-processed material, because even metal- 
poor stars switch to CNO-burning early on. 

Aims. We address the question of whether the CNO surface enhanced self-enrichment in the first few generations of stars could impact 
their mass-loss properties. 

1 Methods. We employ Monte Carlo simulations to establish the local line-force and solve for the momentum equation of the stellar 

(— | _ outflow, testing whether an outflow can actually be established by assessing the net acceleration at the sonic point of the flow. Stellar 

evolution models of rotating metal-poor stars are used to specify the surface chemical composition, focussing on the phases of early 
I ■ enrichment. 

Results. We find that the mass-loss rates of CNO enhanced metal-poor stars are higher than those of non-enriched stars, but they are 
much lower than those rates where the CNO abundance is included in the total abundance Z. Metal-poor stars hotter than ~50000 K, 
in the metallicity range investigated here (with an initial metallicity Z < 10~ 4 ) are found to have no wind, as the high-ionization 
species of the CNO elements have too few strong lines to drive an outflow. We present a heuristic formula that provides mass-loss 
estimates for CNO-dominated winds in relation to scaled-down solar abundances. 

Conclusions. CNO-enriched or not, the winds of metal-poor stars are generally found to be too weak to contribute significantly to the 
overall mass loss. Population III supernovae are thus expected to be responsible for the bulk of the early nucleo-synthetic enrichment, 
unless additional mass-loss mechanisms such as j] Carinae type eruptions or steady mass loss close to the Eddington/Omega limit is 
' important. 

' Key words. -stars:early-type-stars:mass-loss-stars:winds,outflows-stars:evolution 

in; 
Ov 

O. 1. Introduction (Haehnelt et al. 2001; Wyithe & Loeb 2003, e.g.) at redshifts 

above six (Becker et al. 2001). Therefore, it is particularly in- 
■ In the present-day Universe, massive stars lose a substantial frac- teres ting to ponder on the question of how much mass these first 
L; j tion of their initial mass through stellar winds, which they expe- generations of stars might have lost through stellar winds prior 
rience from the moment they form to the moment they end their to explosion 
life as a supernova (SN). These winds are driven by radiation 
J-j j pressure on spectral lines, and their strengths are thus expected The existence of VMS with masses up to 300M o has recently 
to be a function of chemical composition; the lower the metal been claimed for the very luminous stars in the LMC cluster 
content, the weaker the winds. The dependence of mass loss on R136 (Crowther et al. 2010; Bestenlehner et al. 201 1). At high 
metallicity, M(Z), has been studied in detail, both observation- metallicity, mass loss is thought to be strong enough to shed most 
ally (Puis et al. 1996, 1998; Mokiem et al. 2007) and theoreti- of the initial mass, such that the objects may retain only ~3OM , 
cally (Vink et al. 2001; Kudritzki 2002; Krticka & Kubat 2006; and are expected to die as a SN Ic (Heger et al. 2003). However, 
Vink & de Koter 2005; Grafener & Hamann 2008). For a recent as mass loss depends on Z, below a certain threshold metallic- 
review, see e.g. Puis et al. (2008). ity, VMS would be expected to retain most of their initial mass. 

Studies of star formation in the early Universe suggest that If this is indeed the case, some, or maybe even a large number 

the first stars (Pop III) were typically more massive than present- of VMS might be expected to die as pair-instability supernovae 

day objects, with ~10 2 M Q for zero metallicity (Z) objects (PISNe) in the early Universe. The stellar yields of such Pop III 

(Abel et al. 2002; Bromm et al. 2002; Nakamura & Umemura PISNe have been computed by Heger & Woosley (2002). These 

2002; Bromm & Larson 2004; Loeb et al. 2008). Very mas- yields are significantly different from those of "normal" core- 

sive stars (VMS) of zero metallicity are considered to be collapse SNII. PISNe produce orders of magnitude more iron 

the prime candidates for the re-ionization of the universe (up to tens of solar masses) and have a conspicuous yield pat- 
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tern with abundances of even Z elements consistently higher than 
those of oddZ elements (the so-called odd-even effect). 

Observations of carbon-rich chemically extremely metal- 
poor (CEMP) stars have been accumulated over the last two 
decades, but the specific PISN signature remains as yet elu- 
sive (Umeda & Nomoto 2005). It appears that there is a piece 
of the puzzle missing in our understanding of the evolution of 
the first stars. One possibility is that Pop III stars were born with 
masses either higher or lower than the range expected for PISNe 
(140-260M o ), alternatively "something else" must happen dur- 
ing their evolution, which would avoid the expected fate as a 
PISN. This could well be mass loss. 

Mass loss at low Z may also have interesting implications for 
globular-cluster research, as globular clusters display a strong 
abundance scatter amongst many of the light elements (Li, C, 
N, O, Na, Mg, Al). One explanation for these peculiar abun- 
dance patterns involves a first generation of massive stars that 
polluted the cluster before the second generation of stars formed 
(Decressin et al. 2007). Finally, if the time-integrated loss of 
mass is very small, or even negligible, the Pop III stars may 
have left intermediate mass black holes (of masses ~10 2 M ), 
which could have provided the building blocks of the super- 
massive black holes detected in the centers of galaxies today 
(Kawakatu et al. 2005). 

It is clear that there are many reasons to study M(Z). 
Kudritzki (2002) and Krticka & Kubat (2006) found that mas- 
sive zero-metallicity objects may drive only very weak winds 
or no winds at all. Evolutionary predictions for very metal poor 
stars, especially if they are rapidly rotating, show the surfacing 
of sometimes large quantities of primary carbon, nitrogen and 
oxygen produced in situ (e.g Yoon & Langer 2005; Meynet et al. 
2006; Yoon et al. 2006; Hirschi 2007). One may wonder whether 
the surfacing of such material may cause (dramatic) changes in 
the mass-loss behavior of these objects. This question is the topic 
of this study. We present tailored mass-loss predictions for evo- 
lutionary stages of metal-poor stars in which substantial amounts 
of CNO processed material is expected to have surfaced. These 
predictions thus differ from the theoretical studies mentioned in 
the first paragraph, in which scaled solar metallicities are used 
for O-stars. 

The first investigation in this direction was performed 
by Vink & de Koter (2005) who studied the winds of CNO- 
enhanced atmospheres at low iron metallicities (down to Z/Z 
= 10~ 5 ). The enhanced wind driving by CNO-elements discov- 
ered was the reason to study whether primary carbon or nitrogen 
in the atmospheres of extremely metal-poor massive stars might 
boost the mass-loss rates, and potentially prevent the occurrence 
of pair-instability supernovae (Vink 2006). Such complete stel- 
lar disruptions are hypothesized to occur for initial masses in 
the range M = 140 - 26OM if mass loss is negligibly small 
(Heger et al. 2003; Langer et al. 2007). Krticka & Kubat (2009) 
investigated the winds of enriched metal-free stars, describing 
the physics of line-driven winds following Castor et al. (1975, 
hereafter CAK). One interesting aspect of their work is that dur- 
ing the initial hot phase (>50000K) of the evolution of metal- 
poor stars with initial mass above ~30M Q also CNO-enriched 
stars fail to produce a stellar wind. Overall, it was concluded that 
the surfacing of primary CNO does not produce stellar winds 
strong enough to have a major impact on the evolution of metal- 
poor massive stars. 

Our predictions are based on a newly developed 
self-consistent dynamical treatment of stellar winds by 
Miiller & Vink (2008) and Muijres et al. (2012) based on a 
description of the line force that results from a Monte Carlo 



simulation. In Sect. 2 we introduce this method and discuss 
the evolutionary tracks on which the adopted chemical com- 
positions are based. Section 3 presents the results. These are 
discussed in Sect. 4, as well as compared to the results of 
Vink et al. (2001) - often applied in evolutionary calculations - 
and Krticka & Kubat (2009). We end with our conclusions. 



2. Method 

In this section, we present our method for predicting mass loss 
rates for massive stars at very low metallicity. Distinctive for 
these computations is that a tailored chemical composition of 
the surface layers is used to compute the radiative line driv- 
ing. The adopted surface abundance patterns are motivated using 
evolutionary tracks for rotating metal-poor stars. In sections 2.1 
and 2.2 we present and discuss our wind models. Consulted evo- 
lutionary tracks are discussed in Sect. 2.3. The adopted grid is 
presented in Sect. 2.4. 

2.1. Wind models 

To determine mass-loss rates, we use an iterative method, as 
developed by Miiller & Vink (2008) and Muijres et al. (2012), 
in which isa-wind model atmospheres (de Koter et al. 1993) are 
used to describe the model atmosphere and mc-wind models 
(de Koter et al. 1997; Vinketal. 1999) to treat the wind dy- 
namics. The latter uses a Monte Carlo approach based on the 
method developed by Abbott & Lucy (1985). The advantageous 
aspects of this approach are that multiple photon scatterings are 
taken into account, that excitation/ionisation changes of the gas 
are treated self-consistently, and that one can quite easily dis- 
sect the relative contributions of individual species to the line 
force. Essential aspects of the wind dynamics are summarized in 
Sect. 2.2. First, we briefly summarize the main assumptions of 
the model atmospheres. For further details we refer the reader to 
the above references. 

We assume the wind to be homogeneous, spherically sym- 
metric and stationary. Our model atmospheres extend from the 
base of the photosphere to approximately 20 stellar radii. For hy- 
drogen, helium, carbon, nitrogen, oxygen, and silicon, we solve 
the ionization state and the occupation numbers of the levels in 
non-local thermodynamical equilibrium (NLTE). The other ele- 
ments are treated using a modified nebular approximation. The 
radiative transfer in lines is calculated using the Sobolev method 
(Sobolev 1960). 

To compute the line force we make a selection of the 
strongest lines from the line-list in Kurucz & Bell (1995). This 
line list is essentially complete for the ions that play an important 
role in driving the winds studied here. Transitions of the CNO el- 
ements of ionization species C v, N vi and O vi or higher are not 
considered. These ions with one or few bound electrons have 
only few lines in the part of the spectrum that coincides with the 
flux maximum of the stellar radiation. As ionization of all these 
species depends on the radiation field in the Hen continuum 
above 54.4 eV they would remain trace species in our models. 
Shocks in the outflows of massive stars (see Kudritzki & Puis 
(2000) for a discussion/review) may heat a fraction of the gas to 
higher temperatures and produce non-thermal emission at soft 
X-ray wavelength. Such effects may increase the importance of 
transitions of these highly ionized species, although it is unlikely 
that an associated increase in line force is significant. We there- 
fore do not expect that the lack of these lines has an important 
impact on our mass loss predictions. Finally, we note that we 
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miss part of the transitions of N vi. However, for similar reasons 
we do not think that this affects our results in any notable way. 

2.2. Wind dynamics 

The main challenge faced in determining the structure of a line 
driven stellar wind is the strong, but delicate, interaction be- 
tween radiation, line force, and atmospheric structure: The ra- 
diation field and the density (or velocity) structure set the exci- 
tation/ionisation state of the gas. The line force results from the 
structure and the state of the gas, while, in turn, the line force sets 
the structure, therefore indirectly also the radiation field. We try 
to tackle this problem through an iterative process. By assum- 
ing a velocity law and a mass-loss rate we can determine the 
excitation/ionization structure. Next, the line force is calculated. 
Based on a fit to the line force the velocity law is re-determined. 
A model is considered converged once the fit parameters of the 
line force fit function and the mass-loss rate do not change sig- 
nificantly from one iteration step to the next. 

The procedure followed is based on the method described by 
Miiller & Vink (2008) and Muijres et al. (2012), that improve on 
the treatment developed by de Koter et al. (1997) and Vink et al. 
(1999), and applied by Vink et al. (e.g. 2000, 2001) Actually, the 
Miiller & Vink and Muijres et al. studies present two approaches 
to deal with the equation of motion. In one solution, referred 
to as best-/? method, the velocity law is assumed to be a yS-law 
(Lamers & Rogerson 1978). A fit function to the line force is 
constructed that is used to find the most appropriate terminal ve- 
locity y M and measure for the rate of acceleration in the lower 
part of the wind, expressed by p. In the other solution the fit 
function to the line force is used to numerically solve the equa- 
tion of motion, referred to as the hydrodynamical method. Both 
methods are compared in Muijres et al. (2012), where it is found 
that the derived mass-loss rates do not differ more than about 
25 percent and that the terminal wind velocities agree within 20 
percent. 

The hydrodynamical method is dynamically consistent 
but is sensitive to the quality of the fit to the line force at the 
base of the wind, and therefore less stable. In the range of 
parameter space investigated in this article, we anticipate this 
to be more problematic than for the main sequence O-stars 
studied by Muijres et al. (2012). For this reason we apply the 
best ^-method from Miiller & Vink (2008). 

The best /3-method is, strictly speaking, not dynamically con- 
sistent. This warrants caution. Muijres et al. (2012) discuss this 
issue in detail. They identify the solution to be physically cor- 
rect only if the effective gravity gr e ff - i.e. the Newtonian gravity 
corrected for forces due to continuum radiation pressure - ap- 
proximately balances the line force gu ne at the location of the 
sonic point r s , i.e. 



01ineOs) - g e «(r s ) 



(1) 



The sonic point corresponds to the radial distance where the 
flow velocity reaches the local sound speed. We compare the 
line force as simulated by the Monte Carlo code to the effective 
gravity and conclude that if the latter differs from the line force 
by more than 20 percent the wind can not be driven by line radi- 
ation pressure alone. Phrased differently: in this case our model 
does not drive a wind. 

The condition Eq. (1) is more sophisticated than the one ap- 
plied by Krticka & Kubat (2009). In their approach they apply a 
test to establish whether a wind is driven that does not involve 




5.6 5.4 5.2 5 4.8 4.6 4.4 4.2 4 3.8 
Log T eff 

Fig. 1. Evolutionary tracks for metal-poor stars. The green track 
is for Mi n i t = 40 M Q , Z = 10~ 5 and for a rotational velocity 
that is v m t = 555kms _I (Yoonet al. 2006). This model evolves 
homogeneously. The red track is for M; n i t = 85 M , Z = 10 -8 
and v rot = 800 km s (Hirschi 2007). The blue pluses denote 
(L,7eff) combinations for which mass-loss rates have been com- 
puted using abundance patterns that are typical for the surface 
composition of the evolutionary tracks. 



solving the equation of motion, and requires specification of the 
wind mass-loss rate. 

2.3. Evolutionary tracks 

Evolutionary tracks for rotating metal-poor massive stars have 
been computed by Yoon & Langer (2005), Meynet et al. (2006), 
Yoon et al. (2006) and Hirschi (2007). The set of models by 
Yoon & Langer are for initially 12 to 60 M stars. They span the 
metallicity range Z = 10- 5 to2xl0~ 3 and cover initial equatorial 
rotation velocities between zero and 80 percent of the Keplerian 
value. In these tracks, transport of angular momentum is due 
to Eddington-Sweet circulations, shear instability, Goldreich- 
Schubert-Fricke instability, and magnetic torques, and is approx- 
imated by diffusion. In most cases, magnetic torques dominate, 
that are described in the context of the Spruit-Tayler dynamo 
(Spruit 2002), as explained by Heger et al. (2005). The tracks by 
Hirschi (2007) are for initial masses between 20 and 85 M and 
cover metal fractions Z = 10~ 8 to 0.02 and rotational velocities 
from zero to 800 kms~'. Instabilities accounted for by Hirschi 
are meridional circulation (treated as an advective process), and 
secular and dynamical shear instabilities (treated as diffusion). 

Typical examples of the two sets of tracks are shown 
in Fig. 1. Interestingly, the rapidly rotating tracks of 
Yoon & Langer (2005) and Yoon et al. (2006) evolve more or 
less along the zero age main sequence (ZAMS) toward higher 
luminosities. Hence, these stars remain hot and compact. The 
tracks of Hirschi (2007) evolve away from the ZAMS, into the 
cool range of the Hertzsprung-Russell diagram. This striking dif- 
ference between the evolutionary tracks of the two groups is 
probably due to the inclusion of magnetic field in the Yoon et 
al. models. The different implementations of transport processes 
between the models might also have an impact and should be 
investigated properly but this is beyond the scope of this paper. 
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Different prescriptions are used to account for mass loss. 
Yoon & Langer (2005) and Yoon et al. (2006) use results by 
Kudritzki et al. (1989), assuming a metallicity dependence M ~ 
Z 69 from Vinket al. (2001). For Wolf-Rayet phases they use 
results by Hamann et al. (1995), reduced by a factor 10, and a 
metallicity dependence as proposed by Vink & de Koter (2005). 
An ad hoc recipe for the effect of the increased CNO abundances 
on the mass loss is used. Hirschi (2007) adopt the mass-loss 
recipes from Nugis & Lamers (2000) for Wolf Rayet phases. For 
O- and B-star phases, they use Vink et al. (2001). For stars cooler 
than 12 500K they use de Jager et al. (1988). An artifact of the 
recipe of de Jager et al. is that for luminous but cool stars (below 
~ lOkK - 20 kK, depending on luminosity), outside of the range 
of validity of the fitting function, a 20-term Chebychev polyno- 
mial, unrealistically large mass-loss rates are produced (of order 
10~ 2 M yr~'). This rapidly strips stars of their outer envelope, 
causing the tracks of Hirschi (2007) to become strongly enriched 
in this part of the HRD. A feedback occurs, as the metallicity as- 
sociated with the enriched surface causes a high mass loss once 
they evolve back to higher temperatures. This implies that the 
surface enrichment in the blue- ward tracks of Hirschi is strongly 
dependent on this (extreme) mass-loss history. 



2.4. Grid 

We limit the grid of models to the temperature range from 20 kK 
to 50 kK and the luminosity interval from 10 5 5 to 10 6 25 L . The 
reason why we did not try hotter models - although predicted 
by Yoon et al. (2006) - is explained below. Essentially, for such 
high temperatures in combination with a low iron group metal 
content we do not find sufficient line driving to sustain a stellar 
wind. We encounter convergence problems of the model atmo- 
spheres for temperatures close to 10 000 K because convection 
is not accounted for. To avoid convection becoming an issue, 
we adopt 20 kK as the low end of the investigated temperature 
range. 

For the surface abundance of carbon, nitrogen and oxygen 
combined, Zcno, as well as the relative abundances of these 
three elements we assumed values that are typical - but, deliber- 
ately, not tailored - for the predictions of Yoon et al. (2006) and 
Hirschi (2007). Having said this, the abundance choices best fit 
the 85 M e , Z = 10~ 8 and v mt = 800km s -1 track for the case of 
Hirschi and the 40 M @ , Z = 10 5 and v tot = 555 km s -1 track for 
the Yoon et al. case. These are the tracks that are shown in Fig. 1 . 
For the Hirschi case we selected two compositions, one that is 
characteristic for the start of the surfacing of CNO processed ma- 
terial by rotational mixing (termed Case 1) and one that is typical 
for the end of this mixing phase (termed Case 2), but before the 
phase of extreme mass loss as discussed in Sect. 2.3. The stars 
computed by Yoon et al. evolve on a chemically homogeneous 
track and are helium stars at the point they are enriched with hy- 
drogen burning products. We adopt a relative CNO abundance 
that is typical for the start of enrichment (termed Case 4). To sep- 
arate the effects of the H/He ratio from CNO abundance effects 
we introduce an intermediate Case 3 in which the H/He ratio was 
assumed to be equal to the H/He ratio from the Hirschi track. An 
overview of these abundance patterns is provided in Table 1 . 

For all these cases we adopt three values of Zcno, i-e- 0.002, 
0.02 and 0.04. For abundances of other elements, we adopt the 
solar values from Anders & Grevesse (1989) and scale them 
down to the metallicities as listed in Table 1 . We note that these 
other elements do not contribute to the wind driving. 



Table 1. Adopted abundance patterns: X and Y values are always 
consistent with the adopted Zcno, but are given here for the case 
Z CN0 = 0.04. 



Case 


Track type 


X 


Y 


C 


N 


o 


Ziron 


1 


Hirschi 


0.56 


0.40 


42% 


38 % 


20% 


1 x 10~ 8 


2 


Hirschi 


0.56 


0.40 


55 % 


25 % 


20% 


1 x 10~ 8 


3 


Yoon et al. 


0.56 


0.40 


57 % 


10% 


33 % 


1 x 10~ 5 


4 


Yoon et al. 


0.05 


0.91 


57 % 


10 % 


33 % 


1 x 10~ 5 



3. Results 

3. 1 . Mass-loss behavior 

The overall trends in the behavior of the predicted mass-loss 
rates are as follows: For iron-group metallicities below ~ 10~ 4 
(in line with Vinket al. (2001)) and no primary enrichment of 
the atmosphere with carbon, nitrogen and oxygen we find that 
we can not drive an outflow. Only for Zcno ~ 10 3 or higher, the 
CNO elements can drive a stellar wind - if circumstances are fa- 
vorable. For temperatures of 50 kK or higher the CNO elements 
are highly ionized. These ionic species - such as C v, N v and O v 
- have relatively simple atomic level configurations, resulting in 
only few lines near the flux maximum of the emergent spectrum. 
The cumulative radiation pressure of these lines, supplemented 
with a contribution by hydrogen and helium lines that is dimin- 
ishing for increasing temperature, is insufficient to power a wind. 
Therefore, we do not predict stars hotter than ~50 kK to have a 
line-driven outflow. 

Table 3 presents the M, Voo, and /3 predictions for our grid. 
For these models, characterized by a line driving that is on the 
brink of being sufficient to drive a wind, numerical uncertain- 
ties as well as the accuracy of the physical solution requirement 
Eq. (1) may play a role in whether or not a solution is found. 
Typical uncertainties in the (iterative) numerical method itself 
are 0. 1-0.2 dex in the mass-loss rate and ~30 % in the terminal 
velocity. To avoid over-interpretation of individual models, we 
opt to present the results in terms of their overall characteris- 
tics. Moreover, to provide a frame of reference we compare the 
computed mass-loss rates with predictions based on the recipe 
of Vink et al. (2001). These reference rates are referred to as My 
and are also provided in Table 3. 

The Vink et al. recipe has been derived for a chemical com- 
position of all elements apart from hydrogen and helium (i.e. the 
metallicity in the astrophysical context) that is scaled to the solar 
composition as given by Anders & Grevesse (1989). We com- 
pare our mass loss M with My by identifying Zcno of the mod- 
els presented here as the metallicity Z as intended by Vink et al.. 
Clearly, as a general trend we may expect that the Vink et al. 
recipe predicts larger mass-loss rates compared to the values 
computed here, as iron-group elements are more efficient in driv- 
ing the wind through the sonic point than are lines of CNO only. 

The My-recipe is valid for stars that are not too close to 
their Eddington limit. In terms of the ratio of the radiation pres- 
sure on free electrons to the Newtonian gravity, i.e. F e = g e /gN, 
all My predictions are for T e ^ 0.4. An extension of these pre- 
dictions for stars close to their Eddington limit by Vink et al. 
(2011) shows that for F e S> 0.7 the standard My-recipe remains 
unaltered. The largest value of T e in the computations presented 
in this study are T e ~ 0.5 for models in which the luminosity 
log L/Lq = 6.25. We conclude that effects of approaching the 
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Eddington limit do not complicate the comparison of our results 
to the My predictions. 

The outcome of the comparison is shown in Fig. 2, for the 
three different values of Zcno, together with a power-law fit. The 
coefficients of these fits are listed in Table 2. In all cases the 
behavior is approximately 



M(Z CN o) = flM"(Z = ZcNo) 



—mV(Z = Z c 
300 v 



NO 



(2) 



where the mass-loss rates are in units of M yr _I . Note 
that the metallicity Z here concerns the CNO abundance, 
which is disjunct from the mass-loss Z relations provided by 
Vink & de Koter (2005) for WR stars where Z referred to the 
Fe abundance. We point out that the relation is only valid for 
the temperature range 20 kK to ~50kK, and as long as the 
iron-group mass fraction is extremely low, i.e. below a few 
times 1O~ 4 Z , in rough agreement with the mass-loss plateau of 
1O~ 3 Z found by (Vink & de Koter 2005) for late-type WC stars. 
The error in the relation is typically a factor of three. 

Let us appraise some aspects of Eq. (2). For a Vink et al. 
mass loss of 10~ 6 M o yr _1 , our models imply a mass loss that is 
0.3 dex lower for Zcno = 0.002; 0.8 dex lower for Zcno = 0.02, 
and 1.0 dex lower for Zcno = 0.04. So, for progressively higher 
CNO enrichments the discrepancy with My increases. This re- 
flects the fact that the driving power of the CNO lines is topped 
once the relevant driving lines become saturated. An increase of 
the iron-group metal content to high metallicities of 0.02-0.04 
does not suffer from this effect as most of the iron lines that drive 
the wind are relatively weak. Consequently, the mass-loss rate 
of a wind in which the gas composition scales to the solar abun- 
dance pattern continues to rise with increasing metallicity. The 
fact that the main CNO driving lines are saturated in the Zcno 
range explored here implies that the absolute mass-loss rates of 
metal-poor stars enriched with primary CNO must be essentially 
independent of Zcno- Indeed, this is what we find. 

3.2. Hydrogen rich versus helium rich stars 

We find that overall it is more challenging for helium rich stars to 
drive a wind, if all other stellar parameters (notably stellar mass) 
remain unaltered. We identify two reasons for this. First, he- 
lium does not provide as many free electrons as hydrogen does. 
Therefore the contribution of the free electrons to the radiative 
force is less. Second, the contribution of the helium lines to the 
radiative force is at maximum of the order of the contribution 
of the hydrogen lines but in many cases less. The latter effect is 
most relevant for the 20 kK models, where the neutral hydrogen 
fraction (although small in an absolute sense) is the highest. 



Table 2. Fit parameters, expressing our predicted mass-loss rates 
for winds driven by carbon, oxygen and nitrogen lines, in terms 
of the mass-loss rates My as predicted by Vinketal. (2001): 
M = a x My. The coefficients are almost independent of Z, and 
can be summarized according to Eq. (2). 
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-2.51 ±0.31 
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± 0.062 
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± 0.095 
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-2.33 ± 043 


0.660 


± 0.070 
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Fig. 2. Derived mass loss rates versus predictions by Vinket al. 
(2001) (denoted My) for Z = 0.04 (top panel), Z = 0.02 (middle 
panel) and Z = 0.002 (bottom panel). The coefficients of the 
power law fits are given in Table 2. 



The mass loss is found to be independent of H/He ratio. This 
could, however, only be assessed for the case Zcno = 0.002. 
We note that although Cases 3 and 4 are representative for 
the homogeneously evolving tracks of Yoon et al. (2006), they 
need not be representative for strongly enriched modestly rotat- 
ing WNL (that may have Zn > 0.15) or WC stars (that may 
have Zc > 0.50). Models representative for such 'classical' type 
of Wolf-Rayet stars have been computed by Vink & de Koter 
(2005). 

3.3. Terminal velocity behavior 

Figure 3 shows the ratio of the terminal velocity to the effective 
surface escape velocity as a function of Zcno, for all converged 
models in our grid. As for the mass-loss properties we also fo- 



6 



L. Muijres et al.: Mass-loss predictions for evolved very metal-poor massive stars 




-CNO 



Fig. 3. Predicted terminal velocities Voo for all models normal- 
ized to the effective escape velocity at the stellar surface u esc . 
Symbols and color codings indicate temperature; black symbols 
represent the H-poor models. Including the outliers towards high 
Vca/v esc , the average of the H-rich models increases with Zcno, 
from 2.6 to 3.9 to 4.7 for Z CN o = 0.002, 0.02 and 0.04, respec- 
tively. 



cus on the overall trends in the terminal flow characteristics. In 
line-driven winds the nature of the driving lines (both in terms 
of chemical composition and ionization) is reflected in the ra- 
tio Uoo/fesc (see e.g. Castor et al. 1975; Kudritzki & Puis 2000), 
therefore we discuss our results in these terms. 

For the H-rich stars in our grid (Cases 1 through 3) the trend 
is that Voo/Vesc is increasing with Zcno- To be quantitative, the 
ratio increases from 2.7 + 0.7 to 3.9 + 1.1 to 4.7 + 2.4 for Zcno 
= 0.002, 0.02 and 0.04, respectively. We note that the scatter in 
iWi>esc is rather large for the highest metallicity, and that for this 
case the models of 20 000 K only yield a mean u^/v^ that is 
substantially lower, i.e. 3.5 + 1.2. 

The acceleration of the flow in the supersonic part of the 
wind is - as in galactic O- and early B-stars - controlled by 
strong lines of CNO (Vink et al. 1999; Puis et al. 2000). In galac- 
tic OB-stars the observed ratio of terminal velocity to escape 
velocity is ~2.6 for stars hotter than 22 000 K and about 1.3 
for cooler objects (Lamers et al. 1995). Predictions of this ra- 
tio for normal O-stars by Muijres et al. (2012) for the best-/? 
method (see Sect. 2.2) yield values that are typically 30% 
higher than observed. Correcting for this, we still find higher 
outflows speeds. The reasons why the CNO enriched metal-poor 
stars produce higher terminal velocities than Galactic stars in- 
volves a combination of a lower mass-loss rate and a more ef- 
ficient driving in the supersonic part of the outflow, due to the 
high abundances of CNO. Relative changes in the abundances 
of carbon, nitrogen and oxygen (as implied by Cases 1, 2 and 3) 
have a rather modest effect on Vca, given the uncertainty in the 
predictions of the terminal velocity (~30%), this is not signifi- 
cant. 

The converged He-rich models (Case 4) at 30 000 K tend to 
have higher u M /i> esc ratios than do comparable H-rich models. In 
the He-rich case the continuum flux in the region where the main 
driving lines are located is larger, therefore these lines are more 
efficient, leading to larger u M . 



4. Discussion 

4. 1. Comparing to other studies 

Mass loss from luminous metal-free O-type stars enriched with 
primary CNO has been studied by Krticka & Kubat (2009), 
using stellar parameters from the evolutionary tracks for first 
stars, formed from metal-free gas clouds in the early Universe, 
by Marigo et al. (2001, 2003). A quantitative comparison of 
our results for metal-poor stars with those of Krticka & Kubat 
(2009) cannot be made because of the large differences in stel- 
lar parameters and surface Zcno abundances. Yet, we provide a 
rough comparison for a limited number of models. Accounting 
for the uncertainties introduced in their model grid spacing, 
Krticka & Kubat found a similar T e ff boundary for which stars 
in similar luminosity ranges can drive winds: stars must be 
cooler than ~60-50 kK. Accepting a ~ 1 0% difference in the stel- 
lar mass, we find similar mass-loss rates for our model with 
Zcno = 0.02, T e g= 20 kK and log L/L G = 5.9 (see their model 
M500-3): the difference is ~0.2 dex. However, we obtain a termi- 
nal velocity that is 20-40 percent lower, depending on the exact 
choice of chemical patterns. 

We note that there are both similarities and differences in 
our treatment of the stellar atmosphere and wind and those by 
Krticka & Kubat. In both methods the state of the gas is treated 
in NLTE, whilst the Sobolev approximation is used to describe 
the transfer of radiation in spectral lines. Krticka & Kubat, how- 
ever, treat the photosphere - from where the continuum radiation 
originates - and wind separately. They adopt an emergent flux 
taken from H-He model atmospheres, whereas we treat the pho- 
tosphere and wind in a self-consistent, unified way. Moreover, in 
our Monte Carlo treatment we self-consistently account for mul- 
tiple scattering effects, which are ignored in the Krticka & Kubat 
approach. For the type of models discussed here, these differ- 
ences are not likely to cause major differences. 

More important, however, is the difference pointed out in 
Sect. 2.2. Krticka & Kubat apply a test in order to assess whether 
a wind can be sustained that is much simpler than the method 
applied here. The advantage of their criterion is that it allows 
one to a priori determine whether a star can drive a wind. 
However, such a test does not fully capture the feedback be- 
tween line driving and wind density, nor does it fulfill the re- 
quirement of dynamical consistency. For the most luminous ob- 
jects, we can compare models at both 30 kK and 20 kK, again 
for Zcno = 0.02 (their models M999-1 and M999-2). Here we 
find a mass-loss rate that is ~0. 3-0.7 dex lower, while terminal 
velocities are lower by ~20 percent. Here it should be noted that 
Krticka & Kubat adopt a mass of 100 M , while the tracks on 
which we rely suggest ~80M G . 

This necessarily very limited comparison suggests that we 
predict roughly similar mass-loss rates for stars of initially 50- 
75 M Q , but that for more massive stars, in the range of 75- 
1OOM our mass-loss rates are lower by a factor of a few. 

4.2. How to apply M(Z C no) 

Vink et al. (2001) point out that their recipe is valid for a scaled 
solar metallicities as low as 1/30 Z for luminosities less than a 
million solar luminosities, and l/100Z o for more luminous ob- 
jects. Their adopted solar values are from Anders & Grevesse 
(1989), in which Z = 0.019. In exploring the line-driving at 
their lowest metal contents, Vink et al. (2001) already detect an 
increasing relative importance of CNO. Moreover, they point out 
that the occurrence of the bi-stability jump at about 35 000 K at 
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very low metal content due to the recombination of ions of a 
CNO element, i.e. C iv to C ra (Vink 2000). We tested whether 
neglecting the up-turn in mass-loss for stars cooler than 35 000 K 
in the present-day Universe in the My -recipe improved or de- 
graded the correlation presented in Fig. 2 and Table 2. This sig- 
nificantly degraded the correlation. 

To provide a simple recipe to adjust the Vink et al. -recipe for 
models that show primary enrichment, we propose the heuristic 
formula: 

M = M(Zcno) + M Y (Z Mm ), (3) 

where initial is the initial metal content of the star and Zcno 
is the primary CNO abundance. M(Zcno) is given by Eq. (2). 
In the My -recipe, the bi-stability jump should still be applied. 
Equation 3 is valid for T e s < 50 kK. For hotter temperatures 
and abundances below Z ~ a few 10~ 4 (Vinket al. 2001), we 
do not predict stars to lose mass in a line-driven wind. The 
validity of the My formula extends to temperatures as low as 
12500kK. The new Mcno predictions have been computed 
down to 20 000 kK. A modest extrapolation may be applied, al- 
though most certainly not below 15 000kK. The highest Zcno 
abundance in our computations is 0.04. Again, some extrapo- 
lation may be acceptable. However, for chemical compositions 
that are more typical for evolved Wolf-Rayet objects (Zcno ^ 
0.1), we refer to Vink & de Koter (2005). The lowest Zcno abun- 
dance in our computations is 0.002. The minimum Zcno for 
which a wind may be driven is a function of luminosity, prox- 
imity to the Eddington limit, and, to a lesser extent, detailed 
chemical abundance pattern. However, the mass-loss rates as- 
sociated with such low CNO abundances will not impact the 
overall mass-loss and/or angular momentum loss of these stars 
(see Sect. 4.3). This largely removes a practical need to define 
the limit below which metal-poor stars may or may not drive a 
wind. 

We point out that in terms of proximity to the Eddington 
limit, the stars investigated show Eddington factors F e up 
to ~0.5. We therefore do not exclude that metal -poor stars 
very close to their Eddington limit may still drive a wind 
through line radiation pressure, as discussed by Kudritzki 
(2002); Vink & de Koter (2002, 2005); Krticka & Kubat (2006); 
Grafener & Hamann (2008); Vink et al. (201 1). 

4.3. Implications for the evolution of metal-poor stars 

Although the surfacing of CNO-cycle products may make the 
difference in whether or not a metal-poor star can produce a line- 
driven wind or not, the actual mass-loss rate such stars feature is 
only modest. Absolute values of M for the grid of stars presented 
here are in the range 10 6 — 10^ 7 M yr~'. Whilst substantial sur- 
face enrichment may only occur after a significant part of the 
life of metal-poor massive stars, even when simply assuming this 
would happen shortly after reaching the main-sequence the cu- 
mulative mass-loss is only ~1 M for a 30 M star and ~1OM 
for a 80 M star. So, at most a few percent of the initial mass, but 
likely less than that. This is similar to the conclusions reached 
by Vink & de Koter (2005) and Krticka & Kubat (2009). 

Purely by line-driven winds, metal-poor and initially metal- 
free stars may inject at total of the order of 10~ 3 - 10 _1 M Q of 
CNO processed material into their (primordial) surroundings see 
Hirschi (2007) (table 3). This is substantially less than the en- 
richment with CNO and other metals that is to be expected in 
the (pair-instability) supernova that ends their lives. Still, it may 
be too early to conclude that a shift from a top-heavy initial mass 



function (IMF), predicted in simulations of the first star-forming 
regions (Bromm et al. 1999; Nakamura & Umemura 2002), to a 
Salpeter (1955) IMF, is due to dust grains formed from super- 
novae material. This outcome is based on the fact that a line- 
driven wind is the only mass-loss mechanism a metal-poor or 
first star may suffer. Rotation close to critical, perhaps in com- 
bination with approaching the Eddington limit (Meynet et al. 
2006) and/or pulsations (Baraffe et al. 2001) may boost the 
mass-loss rate of these objects. Mass-loss eruptions in a post- 
main sequence phase, as for instance suffered by the massive 
star j] Carinae (Humphreys & Davidson 1994), might play a role 
in the early Universe as well. So, far the exact nature of these 
eruptions eludes us (see Owocki et al. 2004), and possibly this 
mechanism does not depend on the chemical composition. If so, 
this type of mass-loss may have played a relatively important 
role in the Universe early-on. 

5. Conclusions 

We presented mass-loss calculations of metal-poor massive stars 
that, as a result of rotational mixing, have become enriched at 
their surface with products of primary CNO burning. From pre- 
dictions of the evolution of such stars we adopted several charac- 
teristic CNO enriched abundance patterns, and for these chemi- 
cal patterns we performed tailored computations of the mass-loss 
rates driven by radiation pressure on spectral lines. The formal 
requirement to drive a wind in this study is that the gravity bal- 
ances the radiative forces at the sonic point of the flow. 
The main conclusions we obtain are: 

1. The mass-loss predictions My (Vinketal. 2001) have not 
been determined for massive metal-poor (Z S> l/100Z o ) 
stars. Applying this formula, and equalizing the CNO en- 
riched chemical pattern to Z, would overestimate the actual 
mass-loss rate. 

2. On the basis of our new calculations, we provide a heuristic 
formula to estimate the mass-loss rates of such metal poor 
stars, for the parameter range L <: 1O 5 5 L , T e g k, 15 000 K 
and F e 5= 0.5, and we propose to apply this formula (Eq. 3) 
- expressed as a scaling to the oft-used Vinket al. (2001)- 
recipe - in evolutionary predictions of such objects. 

3. Metal -poor and surface CNO enriched massive stars hotter 
than 50 000 K do not feature a line-driven wind; highly ion- 
ized carbon, nitrogen and oxygen ions provide too few ef- 
ficient driving lines. This result is similar to that obtained 
for metal-free stars by Krticka & Kubat (2009), within the 
uncertainties in the T e g resolution of the two grid computa- 
tions. It implies that homogeneously evolving stars (which 
remain very hot during their entire life) do not have line- 
driven winds in the initial metallicity range investigated here 
(Z < 10~ 4 ). 

4. In the range of Zcno = 0.002-0.04, the mass-loss rate is not 
a strong function of the ratio of C to N to O. 

5. The terminal wind velocities of CNO enriched metal-poor 
stars are higher than those of massive stars in the present- 
day Universe. For Zcno ~ 0.02-0.04 they are 25-50% higher. 
The reasons for this behavior are a combination of the lower 
mass-loss rate and a more efficient driving in the supersonic 
part of the outflow, due to the high abundance of CNO. 

6. The winds of massive very metal -poor stars (Z < 10~ 4 ), 
whether they are CNO enriched or not, are so weak that they 
do not significantly impact the total mass and/or angular mo- 
mentum loss during their evolution. If other mass-loss mech- 
anisms, such as T] Carinae type of mass eruptions, do not oc- 
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cur for such objects, their supernova explosions are expected 
to be responsible for the bulk of the early nucleo-synthetic 
enrichment. 
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Table 3. Adopted model parameters and predicted wind properties. Both the models for which a wind solution is achieved and those 
for which this is not the case are listed. The chemical compositions (column 7) correspond to the coding as given in Table. 1 . Typical 
errors on the theoretical mass-loss rates are 0.1-0.2 dex; the errors on the terminal velocities are ~30%. For models in which the 
abundance pattern is followed by a semi colon the errors are larger. The reference mass-loss rates My are as given by Vink et al. 
(2001), if we employ u M =2.6 u esc . Beware that in order to compute My we have set the value of Zcno equal to the metallicity value 
Z, which is not intended by the My-recipe, in which Z refers to the initial iron abundance. Small changes in the effective escape 
velocity from the stellar surface (column 10) are due to abundance effects on T e . 
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-6.416 


1412 


450 


0.81 








0.002 


-6.299 








- 


- 


6.00 65.0 


50.0 


13.3 


0.04 


-4.840 








- 


- 








0.02 


-5.096 








- 


- 








0.002 


-5.946 








- 


- 




40.0 


20.8 


0.04 


-4.828 








- 


- 








0.02 


-5.084 








- 


- 








0.002 


-5.934 


1 


-6.305 


2500 


863 


0.72 












2 


-6.318 


2400 


863 


0.71 












3 


-6.352 


2752 


863 


0.75 




30.0 


37.0 


0.04 


-5.115 


1 


-6.212 


4740 


661 


1.06 








0.02 


-5.371 


1 


-6.160 


2450 


661 


0.78 








0.002 


-6.221 


2 


-6.215 


1305 


660 


0.64 












3 


-6.230 


1205 


659 


0.63 












4 


-6.301 


2500 


736 


0.76 




20.0 


83.3 


0.04 


-4.854 


1 


-6.132 


1625 


449 


0.80 












2 


-6.130 


1350 


449 


0.71 












3 


-6.084 


905 


449 


0.65 








0.02 


-5.110 


1 


-6.236 


1950 


449 


0.96 












2 


-6.144 


1450 


449 


0.82 








0.002 


-5.960 


1 


-6.172 


1250 


449 


0.76 












2: -6.309 


1620 


449 


0.87 
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Table 3. Continued . 



L. Muijres et al.: Mass-loss predictions for evolved very metal-poor massive stars 



Model Parameters 












Wind properties 








L M 


^eff 


R 


ZcNO 


A7 V 


Abund. 


M 


Vet 


"esc 


P 


log(L ) M B 


kK 


Re 




log(M /yr) 


pattern 


log(A/ Q /yr) kms 1 


km s -1 




6.25 80.0 


50.0 


17.8 


0.04 


-4.410 




- 


- 












0.02 


-4.666 




- 


- 












0.002 


-5.516 




- 


- 








40.0 


27.8 


0.04 


-4.398 




- 


- 












0.02 


-4.653 




- 


- 












0.002 


-5.503 


1 


-5.962 


2310 


71 1 


0.76 












2 


-6.004 


2616 


709 


0.80 












3; 


-6.013 


2721 


709 


0.81 












4 


-6.137 


3400 


840 


0.75 




30.0 


49.3 


0.04 


-4.685 


2 


-5.943 


5400 


554 


1.16 








0.02 


-4.941 


2: 


-5.837 


2274 


552 


0.84 












3 


-5.843 


2000 


552 


0.77 








0.002 


-5.791 


2 


-5.879 


1483 


552 


0.74 












4 


-5.927 


1667 


669 


0.65 




20.0 


111.0 


0.04 


-4.422 


1 


-5.786 


1110 


382 


0.74 












2: 


-5.877 


1567 


382 


0.81 












3: 


-5.761 


1240 


382 


0.81 








0.02 


-4.678 


1 


-5.978 


2000 


382 


0.91 












2 


-5.888 


1615 


381 


0.93 












3: 


-5.929 


1600 


381 


0.83 








0.002 


-5.528 













